Magnetism and Transport for Two-Dimensional Electrons : Absence of Stoner 
Ferromagnetism and Positive In-Plane Magnetoresistance in the Metallic Phase 
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We calculate the interaction kernel K for two-dimensional diffusive electrons. The screening of 
the Coulomb interaction together with the Fermi statistics induces a spin selection rule for electron- 
electron scattering so that in leading order in the inverse conductance only pairs of electrons with 
antiparallel spins do scatter. At low temperature, this results in a larger coherence length for fully 
polarized electrons and thus in a positive in-plane magnetoresistance. An applied in-plane magnetic 
field also induces a norimonotonous behavior of K at finite temperature. Alternatively, the vanishing 
of the scattering in the triplet channel strongly reduces ferromagnetism deep in the metallic regime. 
These effects weaken as the density of charge carriers is reduced. 

PACS numbers : 73.23.-b, 75.75.+a, 71.10.-w 



Over the past five years, many experiments have 
reported a nontrivial connection between transport 
and magnetic properties of dilute, high-mobility two- 
dimensional electron gases (2DEG). Especially intrigu- 
ing is the response of those systems to the application 
of a magnetic field parallel to the confinement plane of 
the charge carriers. The magnetoresistance in response 
to such a field is enormous, leading to the suppression 
of the metallic behavior and the saturation of the resis- 
tance above a field strength Bgat 0] ■ Shubnikov-de Haas 
measurements indicate that this saturation occurs as the 
electrons become fully spin polarized while no re- 
sponse was observed for systems with easy-spin axis ||] 
: the effect is related to the electron spin. Recent ex- 
periments report a significant increase of the magnetic 
susceptibility as the conductance of the sample decreases 
and it has even been reported that iJsat vanishes as 
the conductivity crossovers from a metallic to an insulat- 
ing behavior where, extrapolated to T = 0, an arbitrarily 
weak in-plane field seems sufficient to fully polarize the 
2DEG 1^ . Quite naturally, this raises the issue of the 
presence of a ferromagnetic instability as the system be- 
comes localized. 

Strong connections between magnetic ordering and lo- 
calization have been known since the early days of quan- 
tum mechanics: Pairs of weakly overlapping, localized 
electrons can reduce their energy via virtual (kinetic) ex- 
change processes, which together with the Pauli exclu- 
sion, leads to an effective antiferromagnetic coupling 
For strongly overlapping, itinerant electrons on the other 
hand, the interaction energy can be minimized by align- 
ing the spins and this leads to an effective ferromagnetic 
coupling 1^. A Hartree-Fock treatment predicts then a 
ferromagnetic instability when the gain in exchange en- 
ergy counterbalances the loss in kinetic energy due to the 
promotion of electrons to higher unoccupied energy levels 
which leads to the Stoner criterion [B| 



d''rid''r2Vi.2(C(l)^/3(l)V'a(2)^A;(2)) = A (1) 



Here Vi_2 — ^('"1,2 = In. ~ ^^2!) is the interaction poten- 
tial, ?/'a(l) = V'Q(n.) denotes a one-particle wavefunction 
and the average (. . .) is taken over wavefunctions close to 
the Fermi level. A is the level spacing at the Fermi level. 

The threshold (|^) overestimates the tendency toward 
magnetic ordering : it is a perturbation result in the 
interaction which often predicts a threshold where elec- 
trons repel each other so strongly that a perturbation 
theory in the kinetic energy is more adequate. In most 
cases the ferromagnetic order is then destroyed by anti- 
ferromagnetic superexchange HJl^. In other instances, 
J < A and paramagnetism prevails. Remarkable excep- 
tions are provided by materials with orbital degeneracies 
for which ferromagnetism is favored by the absence of 
any kinetic energy cost for spin alignment. This is the 
essence of Hund's first rule and is also believed to lead 
to ferromagnetism in Fe, Co and Ni pT[ . 

Roughly speaking, itinerant electrons have a tendency 
toward ferromagnetism, while the coupling becomes an- 
tiferromagnetic between localized electrons or in the case 
of a strong local repulsion. Below we will establish a re- 
versed connection between magnetic order and localiza- 
tion for disordered systems where ferromagnetism is sup- 
pressed beyond the Thouless energy Ec = gA [g is the 
conductance) in the delocalized phase while the mech- 
anism behind this suppression weakens as the system 
becomes insulating, resulting in an increasing magnetic 
susceptibility in qualitative agreement with experiments 
and numerical simulations ]l^ . 

For diffusive metals, extending the Hartree-Fock treat- 
ment of the Stoner ferromagnetism seems appropriate : 
the (kinetic) antiferromagnetic superexchange is damped 
by (nonmagnetic) impurity scattering. Simultaneously, 
the presence of disorder induces stronger wavefunction 
correlations, effectively increasing the strength of the fer- 
romagnetic exchange (^ whereas it strongly reduces by a 
factor 0(1/(7) ^ 1 the typical magnitude of off-diagonal 
interaction matrix elements ]l3| , p^ . It thus seems rea- 
sonable to truncate the perturbation expansion in the 
interaction. Doing so at first order predicts an increased 
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probability for weak but finite ground-state spin polar- 
ization due to disorder jisl] , a conclusion that has found 
experimental confirmation in ballistic microcavities . 
When one considers higher order corrections in the in- 
teraction one finds however that in either the clean 
or in the disordered case |Q, Hartree-Fock approaches 
overestimate the magnetic ordering. In this letter we will 
amplify on these latter results by a microscopic calcula- 
tion of the relevant parameters of the perturbative ex- 
pansion developed in [|l8| . We will see that in a diffusive 
metal, and assuming a Spin Rotational Symmetric (SRS) 
Coulomb interaction, electrons with parallel spins do not 
scatter in leading order in 1/g, which leads to a strong 
demagnetizing effect in the metallic diffusive regime. 

Our arguments can be summarized as follows. SRS im- 
plies that interaction-induced transitions occur between 
pairs of fermions with conserved total spin (j = or 1). 
Fermi statistics requires the symmetrization (antisym- 
metrization) of the off-diagonal matrix elements connect- 
ing j = (j — 1) two-particle states. For electrons with 
parallel spins, this prescription results in the vanishing of 
the interaction for distances shorter than the Fermi wave- 
length r <kp^ and this strongly reduces off-diagonal in- 
teractions for triplet-paired electrons. The effect becomes 
stronger as kp is reduced: writing the inverse screen- 
ing length as K2 = a2kF, with a2 = l/{kpa) > 1 (a is 
the Bohr radius), one has, {g^ is the valley degeneracy), 
«2 = V^gv^^rs p9| . rg is the gas parameter, hence a2 
increases as fc^ is reduced. 

Three important consequences are : (i) Extending the 
second order perturbation theory developed in Q we 
get a much larger energy cost for polarizing spins and this 
suppresses the Stoner instability beyond the Thouless en- 
ergy. At lower this energy cost is reduced and the 
magnetic susceptibility increases as g ^ 1. This effect 
is enhanced by a simultaneous increase of the ferromag- 
netic coupling constant so that extrapolating our results 
to g ~ 1 (where however the theory loses its validity) 
we reach a ferromagnetic instability, (ii) The strongly 
reduced scattering for electrons with parallel spins in- 
creases the coherence length at low temperature for po- 
larized electrons. As an in-plane magnetic field is applied 
to a 2DEG, this is accompanied by an increase of weak 
localization corrections and thus by positive magnetore- 
sistance. (iii) In presence of a magnetic field, the Zeeman 
splitting Auj = B results in an energy shift for the inter- 
action kernel which becomes: K{uj; T / i) = K{\uj ± B\), 
where the "+" and "-" signs correspond to an electron 
with excitation energy uj being in the majority (t) and 
minority (J,) spin subband respectively. Accordingly, 
the average K = (if (a;; t) + A'(tj; |))/2 shows a non- 
monotonous behavior at finite excitation energy (or tem- 
perature) UJ, increasing with B for B < uj and decreasing 

K± = miim''^* ± (C/J:')*> - 8 / X{d^r,V,.2 



again for B > lu. 

We now present our calculations. Under requirement of 
SRS and rotational symmetry in Hilbert space a generic 
Hamiltonian for interacting electrons in a weakly disor- 
dered metallic system can be written fl^ , ^ 

H = Ho- JS ■ S + Uf. (2) 

The first term is a one-body, spin independent Hamil- 
tonian with eigenfunctions The second term repre- 
sents the ferromagnetic exchange, J > 0. We have in- 
troduced spin operators Sa = (1/2) t Ca,s'^s,tCa,t and 
S — Sa, where the c's are fermionic operators. The 
third term finally contains the physics beyond the mean- 
field approach, i.e. the fluctuating part of the interaction 

^f= J2^^i':ks4,s''^s,s'C^,s, (3) 

a,/9;7,(5 s,s' 

where the interaction matrix elements are defined as 

uli = / d'nd\2Vi,2ram;msmM'^)- (4) 

We will use the d = 2 screened Coulomb potential 

Vi.2 = {e/2nf J d^qexp{iq- ri^2)/{q + K2{q,io)), (5) 

where screening depends on the energy transfer lu, 
k2{q,uj = 0) = K2 [0 and we wrote fi^2 — ri ~ f2- 

Uf can be rewritten as a sum over four contribu- 
tions : the first one inducing scattering between singlet- 
paired {j = 0) and the other three between triplet-paired 
(j = 1, M = 0, ±1) two-electron states, defining electron- 
electron scattering in singlet and triplet channels respec- 
tively. Accordingly, the interaction matrix elements must 
be symmetrized or antisymmetrized corresponding to the 
singlet ("-I-" sign in the formula below) and the triplet 
("— " sign) channel respectively 

From (^) and (^) , the interaction matrix elements in both 
channels are symmetrically distributed around zero aver- 
age and the interaction kernel is related to the variance 
of these distributions K± = a'^{U^). Uf systematically 
opposes ferromagnetism, as it induces more virtual pro- 
cesses in low-spin sectors [ p^ . Since in the ferromagnetic 
sector perturbative corrections to the energy contain only 
contributions from the triplet channel, this tendency is 
further amplified if ^ K+. We will now show that 
this is indeed the case in a metal. 

After disorder averaging, K+ and AT- are given by 

,,4 [(C(l)V'7(l)'/'o(4)V'; (4)) (V'^ (2)^5 (2)^-/3(3)^1(3)) 
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(7) 



For diffusive systems, 1^ ^ L [l^ is the elastic mean free 
path, L the linear system size), one has up to 0(1/5) 11] 

(C(l)V'7(l)V'a(4)V;(4)) = n(ri,4)/f^^ (8) 

il = L'^ and the diffusion propagator is given by n(r) = 
l/ii-K^iy) J d'^qexp{iq-f)/{Dq^-iLo). D = gl{2-nv) is the 
classical diffusion constant, v is the density of states, the 
integral ranges from to and uj gives the energy 
difference between the states a and 7. When r3_4 < 
the second term between brackets in the integrand of (m) 
is equal to the first one and we can readily conclude that 
for a strong screening a2 ^ 1 the triplet channel scat- 
tering is completely suppressed (this limit corresponds to 
the Hubbard model). For r3^4 > k'p^ , one has 

(C(l)V'7(l)^a(3)V;(4)) (X exp(-(ri,3 + ri,4)/(2?e)), (9) 

and (||) vanishes exponentially for 1^ L. We write 
K± = K ± SK. A straightforward calculation gives K — 
2(A/7r2g)2. To calculate we note that, in the strongly 
diffusive regime the diffusion modes are strongly damped 
for q > Z^^, furthermore one has K2le ^ kpl^ S> 1. This 
justifies an expansion of the Coulomb potential (js]) in in- 
verse powers of K2 ■ Keeping only the first two nonvanish- 
ing orders (which is consistent with the above expression 
for the diffusion propagator and the truncation of (|^) at 
the first order in 1/g) one gets for lo < Ec 

SK = K{1 - 4ia2gy\le/L)) + 0{g-^), (10) 

and for uj > Ec, both K and SK must be multiplied 
by a factor (1 +uj/Ec)^^ Q. K^ is thus suppressed 
by a factor {a2g)^^le/ L <C 1 compared to K-^- which 
reflects the fact that the long-range part of the RPA po- 
tential (1) gives subdominant contributions to K as long 
as K2le ^ 1, and we thus set K^ — 0. Note that for 
d = 3, the discrepancy between K^ and K- is ampli- 
fied by the exponential screening of the interaction while 
balhstic systems have K+/K^ = 0(1) ||2|]. 

Up to 0(1/ g), the average exchange strength (^ reads 

J^{A/2)[Jsr + Jb + Jd] (11) 

JsR, Jb and Jd are short-range (for ri_2 < 1/^2 in (0)), 
long range (ri_2 > 1/^*2) ballistic and diffusive contribu- 
tions respectively 

JsB = l + \og{L/Q/{T:g) 
Jb^ dRJ^{R)/{a2R^) (12) 

Jd - (1 + JB)log(L/g/(^5) 

For kple — g ^ 1, JB{a2) ~ 1 — (1 — 7r^^)/a2. Inter- 
estingly for Q!2 3> 1, one has J > A and a first order 
treatment predicts a ferromagnetic (Stoner) instability. 



We can now calculate the average ground-state energy 
in each spin sector up to the second order in the interac- 
tion. For fixed 5*, the number of majority (-I-) and minor- 
ity (-) spins are n± = n/2±S where n is the total number 
of electrons. We first neglect energy dependences in (|[]|) 
which reduce K± and accordingly truncate the Hamilto- 
nian at the Thouless energy, i.e. we set n — g. The 
number of singlet (triplet) transitions in leading order is 
then N+ = g'^n+n^/A {N^ 9^{n\ +11"^+ n+n_)/4) 
where the prefactors g^ account for the final states. In 
presence of a magnetic field B, the average energy dif- 
ference between the paramagnetic ground-state and the 
ground-state in the ^-sector in second order perturbation 
theory is given by 

As = JSiS + 1) + BS-S^A + g\og{g)SK/{2A)] (13) 

which illustrates the demagnetizing influence of SK. For 
large metallic samples with 5 3> 1 and 1^ <ti L one 
has J < A and the magnetization is bounded by Sq < 
Smax — 0{g/ \og{g)). Electrons are not polarized beyond 
the Thouless energy Ec and hence there is no thermo- 
dynamic magnetization. At the same time this does not 
rule out the probability to flnd few polarized electrons 
in agreement with experiments on quantum dots in the 
Coulomb blockade regime (In that latter case how- 
ever, ballistic contributions that were neglected in (|^) 
have to be considered |2^].) At lower g on the other 
hand, one may get J > A (a2 increases), and since gen- 
erally log{L/le) > log(f/) for d = 2, the polarization does 
not saturate below Ec- For larger energies u! > Ec, be- 
cause of the decay of K, second order corrections become 
marginal, and therefore the polarization becomes ther- 
modynamical. While this latter point should be taken 
with a grain of salt (even nonmagnetic materials may 
have an exchange above the Stoner threshold (0), see dis- 
cussion in the introduction) , this argument demonstrates 
that the magnetic susceptibility increases as g decreases. 

We next calculate the magnetic field Bgat necessary to 
achieve full polarization (ASgat = S^at - B°^^; B°^^ = 
2g{A — J) — J is the first-order saturating field). For 
5 3> 1, Ai?sat is dominantly given by the field it takes to 
polarize the electrons up to the Thouless energy: 

Ai3,at « 4A(1 - 4{a2g)-\lc/L)) \og{g)/TT\ (14) 

Clearly, Bgat decreases when g is reduced in agreement 
with the experiments U,^. 

We finally calculate the electron-electron scattering 
rate for an electron at an excitation energy uj above a 
Fermi sea polarized by an external magnetic field B. 
For an electron in a given spin subband (the + and — 
signs correspond to the majority and minority subband 
respectively) this rate is determined by the scattering 
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with a particle below the Fermi energy in the other spin 
subband with an energetic difference of w ± i?. Conse- 
quently, K becomes spin-dependent, K{uj) K{iLj; t / i 
) = kIo){1 + \u)±B\/Ec)'^^^-^ and for d = 2 one gets an 
average scattering rate F = (F-|_ -f F_)/2 with 

r±{B,uj)^i27T/h)K{uj;U lysiu^) (15) 
= {2TT/h)K{0)iy3{u)il + \iu± B\/E,)-^ (16) 

U3 is the three-particle density of states which we assume 
is spin-independent. For B > Ec, one then has a non- 



monotonous behavior of F(_B,cl') vs. uj which in particu- 
lar develops a local maximum a,t uj = B for B > 2.28Ec- 
This is so because of the nonmonotonous behavior of the 
average interaction kernel = {K{uj;'\)+K{uj; \,))/2) 



dK{u)/duj>0 ; lu < B, 
dK{Lu)/dLu<0 ; lu>B. 



(17) 



At low temperature {T = lu), dephasing is primarily due 
to electron-electron scattering, and for u!,B ^ Ec the 
weak-localization correction to the conductivity read 



Sae/Sao = 1 + {log[L^/2] + log[l/(L^ + B) + l/\u^ B\]}/ log[F(0, u;)n 



(18) 



where is the elastic mean free time. Equations (|l5| - 
[isl ) predict a nonmonotonous behavior of the quantum 
corrections to the conductivity as a function of an in- 
plane magnetic field. Accordingly, the conductivity first 
increases at low field until B — uj and then decreases 
for larger fields. This should dominate the behavior of 
the conductance for low enough temperature and good 
metallic samples. 

Experimentally, the interaction kernel K can be ob- 
tained via measurements of the energy distribution of 
quasiparticles as e.g. in metallic quantum wires ||2^ ]. The 
geometry of those systems differs from the twodimen- 
sional case studied here, resulting in a different energetic 
decay oiK{uj) . However as long as the magnetic field 
induces only a Zceman coupling one expects 



K{B,uj) cx {lu + B)/Ec 



\uj-B\/Ec 



(19) 



Recently, such a behavior has been observed for cop- 
per wires |23|. While the presented theory does not ex- 
plain why experimentally a ^ —3/2 as expected [ p^ , 
( p^ ) nicely reproduce the experimentally obtained non- 
monotonous behavior of K{uj) (see Fig. 5 in [^). Al- 
ternatively, the nonmonotonicity of K{oj) could be due 
to magnetic impurities (which would also explain the 
anomalous behavior K{bj) oc (Ec/oj)"' with a ^ d/2 — 2) 
|^,|^, therefore similar experiments on silver samples 
where a — d/2 — 2 and thus magnetic impurities do not 
play a dominant role |^ would be highly desirable. A 
detailed calculation for the geometry corresponding to 
those experiments will be published elsewhere ||2l[| . 

In summary we have found that spin selection rules 
for electron-electron scattering in diffusive metals induce 
a nontrivial coupling between transport and magnetic 
properties of 2DEG and a nonmonotonic behavior of the 
dephasing rate in presence of a Zeeman coupling. 
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